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Abstract 

Most of the reliability literature on modeling the effect of repairs on systems as- 
sumes the failure rate functions are monotonically increasing. For systems with 
non-monotonic failure rate functions, most models deal with minimal repairs 
(which do not affect the working condition of the system) or replacements (which 
return the working condition to that of a new and identical system). We explore 
a new approach to model repairs of a system with a non-monotonic failure rate 
function; in particular, we consider systems with a bathtub-shaped failure rate 
function. We propose a repair model specified in terms of modifications to the 
virtual age function of the system, while preserving the usual definitions of the 
types of repair (minimal, imperfect and perfect repairs) and distinguishing be- 
tween perfect repair and replacement. In addition, we provide a numerical illus- 
tration of the proposed repair model. 



1 Introduction 

Most engineered systems - defined as an arrangement of components that to- 
gether perform an identified (and predefined) set of functions - are susceptible to 
failures, and require some form of rectification in order to return to a functioning 
condition. Most rectifications have an effect on the probability and number of 
future failures of the system over a given period of time. 

The sequence of numbers of failures of the system in time, i.e. the failure pro- 
cess, is modeled as a stochastic counting process, assuming that there can be at 



most one failure in an infinitesimally small interval of time. When the rectifica- 
tion action following each failure is immediate and instantaneous, this process 
can also be described as the sequence of times to failure (or consecutive system's 
lifetimes). 

A counting process is completely described by its conditional intensity func- 
tion, and therefore, rectifications are usually defined in terms of their effect on the 
conditional intensity function of the failure process. The initial conditional inten- 
sity function is the failure rate function of the original lifetime (time to first failure 
of the system), which is often a continuous function of time, and is classified as 
constant, monotonic increasing or decreasing, or a combination of these. Beyond 
the first failure, the conditional intensity function is altered in accordance with 
the rectifications performed following the first and all consequent failures. 

Not all rectifications have the same effect on the system, and based on their 
effect, they are categorized as either replacements or repairs with varying degrees 
of effectiveness. In some cases, replacements can be viewed as extreme repairs. 

In this article, we suggest an approach to model the effect of rectifications 
(here, repairs) for a system having a non-monotonic (here, bathtub-shaped) fail- 
ure rate function. We define repairs in terms of their effect on the virtual age of 
the system. 

The article is arranged as follows. In Section |2] we discuss the concepts men- 
tioned above in more detail, and provide a brief review of existing models rele- 
vant to our study. In Section |3l we describe the repair model and provide model 
formulation. In Section HI we provide a numerical illustration of the proposed 
model. Finally, in Section |5l we conclude with a discussion of the proposed model 
and some directions for future research. 



2 Background and Definitions 

Let X c {i) denote the conditional intensity function of the failure process denoted 
by {N(t);t > 0}. Then 
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where N(t + dt) — N(t) is the number of failures in the interval (t, t + dt], and 
T t = {N(s); < s < t} is the history of the process before time t. The initial 
conditional intensity (or baseline intensity), denoted by A (t), is the failure rate of 
the time to first failure, which is 
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A failure rate or the corresponding distribution is categorized as: constant failure 
rate (CFR) when r(i) is constant over t; increasing failure rate (IFR) when r(t) is 
increasing in t; decreasing failure rate (DFR) when r(t) is decreasing in t; or some 
combination of these. For instance, the bathtub-shaped failure rate (BFR) function 
which is initially decreasing, then constant and finally increasing: 

nit) : r[(t) < , t<ai 
r(t) = { r 2 (t) : r' 2 (t) = , a x < t < a 2 (1) 
r 3 (t) : r' 3 {t) > , t > a 2 , 



where a\ and a 2 are the change points (points at which the the derivative of the 
failure rate function changes sign) of the BFR function. The BFR function is a 
generalization of the above categories; setting oi = a 2 = (ai = a 2 = oo or ai = 
and a 2 = oo), it becomes an increasing (decreasing or constant) function. Setting 
a± = a<i = a, we get a U-shaped failure rate (UFR) function (|Lai and Xie 2006"]) ; see 
Figure CD 
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(iii) (iv) 
Figure 1: Failure rate functions: (i) IFR; (ii) DFR; (iii) UFR; (iv) BFR. 



The quantity r(t) dt is the approximate probability that the system will fail for 
the first time in (t, t + dt]. The distribution of the time to first failure, denoted by 
Ti, is defined by the failure rate r(t), but subsequent failure times are affected by 
the type of rectification performed after a failure. 

Rectifications are broadly classified into repair and replacement. With replace- 
ments, the system is replaced by a new, identical system upon failure. The condi- 
tion of the system following a replacement is therefore identical to a new system. 
In this case, the failure process is modeled as a renewal process with conditional 
intensity function X c (t) — r(t — 2V(t)) (where T N ^ is the time of the last replace- 
ment (perfect repair)), and the expected number of replacements is given by the 
renewal function (jHokstad 1997[ IHunter 19741 |Aven and Jensen 1998| ). Replace- 
ments apply when the system is beyond repair (or non-repairable) or when re- 
placing the system is more feasible than repairing it. 

Repairs are characterized by their effectiveness, which is often referred to as 
the degree of repair - the degree to which the functioning condition of the system 



is restored following the repair. Based on this degree, repairs are typically cate- 
gorized as minimal, imperfect or perfect repair. 

Minimal repairs have no effect on the failure rate function; i.e., the system im- 
mediately before and after the failure is in the same functioning condition. There- 
fore, when all repairs are minimal, the failure process is a non-stationary Poisson 

process with conditional intensity function X c (t) = r(t). The expected number of 

t 

failures over an interval [0, t) is then given by E[N(t)\ = J A c (s) ds. 

o 

Perfect repairs, for IFR functions, leave the repaired system in an as-good-as- 
new working condition, which implies that a perfect repair is equivalent to a 
replacement and can be modeled as one. This definition works for systems having 
an IFR function, since the system at the start of its lifetime has the lowest failure 
rate, i.e., it is at its best. However, if a system has a failure rate function that 
is initially decreasing (e.g. BFR), this definition does not hold, because a repair 
that leaves the system in an as-good-as-new condition is actually worsening the 
system, since the failure rate of the system at the start of its lifetime is higher than 
its failure rate when it is working at its best. Therefore, we distinguish between 
replacement and perfect repair, and describe a perfect repair as the best form 
of repair (not taking into account improvements /upgrades). In other words, a 
perfect repair is one that restores the functioning condition of the system to its 
condition when it is performing at its best. This point of ideal performance is at 
the start of the system's lifetime for a system with an IFR function, but not for a 
system having an initially decreasing failure rate function. 

Imperfect repair, sometimes referred to as general repair, is any repair that 
leaves the system in a functioning condition that is between the functioning con- 
ditions following minimal and perfect repairs. For systems with a IFR function, 
the definition of an imperfect repair has included the extremes minimal repair 
and replacement (aka perfect repair). Here, imperfect repair includes as its ex- 
tremes minimal and perfect repairs, but not replacements. Therefore, we distin- 
guish between repair and replacement. 

In most settings, the degree of repair is a variable with range [0, 1], where a 
degree of zero corresponds to a minimal repair, a degree of one corresponds to a 
perfect repair, and a degree between these extremes corresponds to an imperfect 
repair. Therefore, the higher the degree of repair, the bigger the improvement in 
the functioning condition of the system. 

Here, we do not consider repairs that can worsen the system or upgrades (or 
improvements) . 

Many repair models have been suggested for systems having IFR functions. 
Some common models are the virtual age models discussed in Kijima (1989[ ), 
Varnosafaderani and Chukova (2012a[ ) and Doyen and Gaudoin (2004[); and the 



intensity reduction models discussed in Lindqvist (1998 ), Varnosafaderani and Chukova (2012b ) 
and Doyen and Gaudoin (2004). Models of repair in the case of BFR functions as- 
sume that rectifications are either minimal or replacements. The virtual age mod- 
els for IFR functions have been applied to BFR functions; see ( Dijoux 2009| ), but 
due to the failure rate being initially decreasing, repairs of degree greater than 
zero actually worsen the product. 

In this article, we propose a new approach to modeling imperfect repairs for 
systems having BFR functions which better suit the definitions of the types of 



repair. The effects of repairs are described in terms of modifications to the virtual 
age function of the system. 



3 Modeling the Effect of Repairs 

Let Tj denote the time of the ith failure (also repair, since repairs are immediate 
and instantaneous), and let 5i denote the degree of the ith repair. Also, let A(t) 
denote the virtual age of the system at time t. 

Based on the virtual age function of the system, we propose the following 
repair model. The virtual age of the system at time t, is given by 



A(t) 
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where A(Ti) is the virtual age of the system at the time of its ith failure. Before the 
first failure, when N(t~) = 0, the virtual age is simply A(t) = t; and immediately 
after oi and before the first failure in the useful life period, the virtual age of the 
system is again A(t) = t. See Figure|2]for an illustration of the virtual age function 
for three failures occurring at times t\, t 2 and t 3 . 

With a failure rate function that is initially decreasing, the point of best perfor- 
mance is not the start of the system's lifetime, but the point at which the failure 
rate function is at its lowest. This point for a BFR function is the first change point 
ax. 



A(t) 
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Figure 2: Virtual age function following imperfect repairs of varying degree. 



Therefore, according to this model, when the virtual age of the system at the 
time of the zth failure is less than the first change point a\, the effect of a repair is 



modeled as an increase in the virtual age of the system, such that, a perfect repair 
results in the virtual age being a\. At a\, the virtual age of the system is set to its 
calendar age, i.e. A(a{) = a\. This extends the useful life period of the system, 
which will decrease the probability of future failures. When the age of the system 
is greater than the first change point a\, then the effect of a repair is a decrease in 
the virtual age of the system, such that, a perfect repair results in the virtual age 
being a±. The point a± is the point of ideal performance, because it is the start of 
the useful life of the system, and the failure rate of the system at this point is at 
its lowest. 

The conditional intensity function of the failure process is given by 

W.N / Ht) , t < 2l 

M ) ~ I Xo(A(t)) , t>T 1 

where Ao(.) is the baseline intensity (or failure rate) function. See Figure |3] for an 
illustration of this function following repairs of varying degree. 
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Figure 3: Conditional intensity function following: four imperfect repairs of varying de- 
gree (dashed line); minimal repairs (solid line). 

The repair model stays true to the definitions of the types of repair. A perfect 
repair is the best form of repair, and should result in the system performing at 
its best, which is in this case at a\. A minimal repair, should by definition leave 
the system in the same condition that it was prior to failure, and here, the virtual 
age does not change following a minimal repair. The effect of an imperfect repair 
should be between those of the minimal and perfect repairs, and effectiveness of 
the repair should increase with its degree. Here, as the degree of repair increases, 



so does the effectiveness of the repair (which is reflected in the decrease in the 
conditional intensity function of the process). 

The assumption for this model is that the useful life period (a 1; a 2 ] of the sys- 
tem is at least as long as the DFR period (0, 01], i.e. a 2 — ai > a±. 

4 Numerical Illustration 

In this section, we provide a simple example that illustrates the proposed repair 
model. 

The baseline intensity function used in this example is 

A + at (ai - t) 131 , t < ax 
X (t) = { X . ai <t<a 2 ( 3 ) 

A + a 2 (k - a 2 ) /32 , t > a 2 , 

where A > 0, /3 2 > 0, j3± > f3 2 , and a\, a 2 > 0. The parameter values are chosen 
to be A = 1, a.\ = 0.6, a 2 = 0.5, /?i = 2.5, and f3 2 = 2.8, and the change points are 
chosen to be a\ = 4 and a 2 = 8. 

Since virtual age models for IFR functions have been frequently examined and 
the effect of repairs in this case is known, we limit our illustration to exploring 
the effect of repairs based on our virtual age model in the DFR phase. To do so, 
we select an arbitrary mission time r, and applying repairs of varying degree in 
the interval [0, ai), we compute the expected number of failures in (0, r]. Here, 
the mission time is chosen to be r = 10. 

The repairs performed are chosen according to the following strategy: the first 
repair in the interval (0, a\\ is imperfect, and all other repairs are minimal. 

Let T\ denote the time of the first failure. The density function of T\ in terms 
of the baseline intensity function is given by 

i 

— f Ao(s) ds 

h(t) = X (t)e o . (4) 

The expected number of failures in the interval [0, r), is then derived as fol- 
lows: 
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where 5\ is the degree of the imperfect repair performed in (0, a x ]. 

Tabulated in Tableware the expected numbers of failures E[N(10)] for degrees 
of repair 5i e {0.1, 0.2, . . . , 1.0}. 

Table 1: Expected number of failures in the interval [0, r) for various degrees of repair 



Si 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


E[N (10)] 


33.78 


27.3 


22.4 


18.81 


16.29 


14.64 


13.63 


13.09 


12.86 


12.79 


12.78 



Note that, according to the repair model, as the degree of repair increases, the 
expected number of failures decreases; also see Figure |U 
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Figure 4: Expected number of failures £"[^(10)] for 5i € [0, 1]. 

5 Conclusion 

In this article, we proposed a new repair model for systems having a BFR func- 
tion. The effect of repairs was modeled as a modification in the virtual age of the 
system following the repairs. 

According to the proposed model (illustrated in Section H]), as the degree of 
any given repair increases (while others remain fixed), the expected number of 
failures decreases, since the reliability of the system is improved. 

Some possible future research directions are deriving virtual age models for 
systems with more than two change points and extension of these models to two 
dimensions. 
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